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Abstract 

Bose-Einstein correlations in pairs of identical charged pions produced in a sample of 4.3 mil- 
lion T? hadronic decays are studied as a function of the three components of the momentum 
difference, transverse ("out" and "side") and longitudinal with respect to the thrust direction 
of the event. A significant difference between the transverse, r^^.^^, and longitudinal, r/, di- 
mensions is observed, indicating that the emitting source of identical pions, as observed in 
the Longitudinally CoMoving System, has an elongated shape. This is observed with a va- 
riety of selection techniques. Specifically, the values of the parameters obtained by fitting 
the extended Goldhaber parametrisation to the correlation function C = qdata^qmc t^^^ 
two-jet events, selected with the Durham algorithm and resolution parameter y^ut = 0.04, are 
r^^^^^ = (0.809 ± 0.009 (stat) t'^ofi (syst)) fm, r, = (0.989 ± 0.011 (stat) t'^omt (syst)) fm and 
ii/it^,i^ = 1.222 ± 0.027 [stat) lo;oi2 {syst). The results are discussed in the context of a recent 
model of Bose-Einstein correlations based on string fragmentation. 
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1 Introduction 



Bose-Einstein correlations (BECs) in pairs of identical bosons, mainly Tr^'^vr^, have been 
widely studied at various energies for hadronic final states produced by different initial states: 
e"'"e~ 1^, ep [§], pp 0, vrp, K^'^p and heavy ion collisions 0. Two-particle BECs have also 
been studied for pairs [Q, for K^K^ [^] and, at LEP2, for pions coming from W~^W~ 



decays |10|. Genuine BECs have also been observed for three charged identical pions [11 



BECs are manifested as enhancements in the production of identical bosons which are close 
to one another in phase space. They can be analysed in terms of the correlation function 

C{pi,P2) = —, r, (1) 

Po(Pi,P2) 

where pi and p2 are the four- momenta of the two bosons, p{pi,P2) is the measured density 
of the two identical bosons and Po{pi,P2) is the two-particle density in the absence of BECs. 
The choice of the reference sample used to determine Po{pi,P2) is crucial for the measurement. 
It should have the same properties as the sample used for p{pi,P2) except for the presence of 
BECs. In this paper we use pairs of particles with charges of opposite sign as the reference 
sample. This sample, however, includes pairs coming from resonance decays and from weakly 
decaying particles, like the K^. In addition, the correlation function obtained with this reference 
sample has to be normalized and suffers, at large four-momentum differences, from long-range 
correlations due to energy, momentum and charge conservation. We therefore also use a large 
sample of Monte Carlo events without the simulation of Bose-Einstein effects in order to obtain 
a correlation function which is self normalized and which has a reduced contamination from 
correlated unlike-charge pairs. This correlation function, called C'{pi,p2), is used in the paper 
to obtain the reference results. 

The information obtained from the shape of the correlation function may be used to infer 
the space-time extent of the particle emitting region. Most analyses have been performed as- 
suming a spherical emitter but several theoretical investigations have recently treated the shape 



of the correlation function in more than one dimension The Lund group, in particular, has 
developed a model for BECs based on a quantum mechanical interpretation of the string area 
fragmentation probability ||13|. One of the main predictions of the model, developed for two-jet 



events, is that, since momentum components longitudinal and transverse with respect to the 
string direction (i.e. event direction) are generated by different mechanisms, the correlation 
length in the longitudinal direction is different from that in the transverse one. In particular, 
in the so called Longitudinally CoMoving System (LCMS), the longitudinal source dimension 
is predicted to be larger than the transverse dimension. 

In this paper we describe an experimental study of transverse and longitudinal BECs per- 
formed with the high statistics sample of hadronic events recorded by the OPAL detector at LEP 
at centre-of-mass energies at and near the Z° resonance. Similar two- and three-dimensional 
studies have been done in e"'"e~ annihilations at centre-of-mass energies of 34 GeV M] and 



91 GeV [|T^. We also present the results of a unidimensional analysis of BECs. 
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2 Detector and Data Selection 



A detailed description of the OPAL detector can be found in refs. []16|, ^7j. This analysis is 
based mainly on the reconstruction of charged particle trajectories and momenta in the central 
tracking chambers and on energy deposits ( "clusters" ) in the electromagnetic calorimeters. All 
tracking systems are located inside a solenoidal magnet which provides a uniform magnetic field 
of 0.435 T along the beam axis The magnet is surrounded by a lead-glass electromagnetic 
calorimeter and a hadron calorimeter of the sampling type. Outside the hadron calorimeter, 
the detector is surrounded by a system of muon chambers. There are similar layers of detectors 
in the barrel (|cos^| < 0.82) and endcap (|cos6'| > 0.81) regions. The central tracking detector 
consists of a silicon micro- vertex detector [|1^, close to the beam pipe, and three drift chamber 
devices: the vertex detector, a large jet chamber, and surrounding z-chambers. The vertex 
chamber is a cylindrical drift chamber covering a range of | cos6'| < 0.95 with a resolution of 
50 /zm in the r0 plane and 700 fim in the z direction. The jet chamber is a cylindrical drift 
chamber with an inner radius of 25 cm, an outer radius of 185 cm, and a length of about 4 m. 
Its spatial resolution is about 135 /im in the r0 plane from drift time information and about 
6 cm in the z direction from charge division. The z-chambers provide a more accurate z mea- 
surement, with a resolution of about 300 /im. In combination, the three drift chambers yield 



a momentum resolution of (Jpjpt ~ ^0.02"^ + (0.0015 ■ Pt)"^ for | cos(^)| < 0.7, where pt is the 
transverse momentum in GeV/c. Electromagnetic energy is measured by lead-glass calorime- 
ters surrounding the solenoid magnet coil. They consist of a barrel and two endcap arrays with 
a total of 11704 lead-glass blocks covering a range of | cos6\ < 0.98. 

A number of selection cuts are applied to the initial data sample, consisting of 4.3-10^ 
hadronic events from Z° decays. In order to be considered, a charged track is required to have 
a minimum of 20 hits in the jet chamber, a minimum transverse momentum of 150 MeV/c and 
a maximum momentum p of 65 GeV/c. Clusters in the electromagnetic calorimeter are used 
in the jet-finding algorithm if their energies exceed 100 MeV in the barrel, or 200 MeV in the 
endcaps. Only events that are well contained in the detector are accepted, by requiring that 
\cos9thrust\ < 0.9, where Othrust is the polar angle of the thrust axis, computed using charged 
tracks and electromagnetic clusters that passed the above cuts. A second set of cuts, specific to 
the BEG analysis, is then applied. Tracks are required to have a total momentum p < 40 GeV/c 
and to come from the interaction vertex. Electron-positron pairs from photon conversions are 
rejected. Events are selected if they contain a minimum number of five charged tracks and if 
reasonably balanced in charge, i.e. if — 'n-chl I {''^th + ^ch) — 0.4, where n^ and n~^ are 
the number of positive and negative charged tracks, respectively. The analysis is performed on 
the inclusive sample of all events passing the above cuts. In order to compare the experimental 
results with the predictions of |]IB[, which are relevant to two-jet events, and to study the 



dependence of the results on the "jettyness" of the event the same analysis is done on samples 
of events defined as "two-jet events" by the Durham jet-finding algorithm [|l^], with various 
values of the resolution parameter jcut- 



^The coordinate system is defined so that z is the coordinate parallel to the e"^ and e~ beams, with positive 
direction along the e~ beam; r is the coordinate normal to the beam axis, (j) is the azimuthal angle and 9 is the 
polar angle with respect to +z. 
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Figure 1: The LCMS frame of reference drawn for a pair of particles where pit and p2t are the 
projections of the two particles momenta onto the plane perpendicular to the thrust axis. 



3 The Longitudinally CoMoving System 



In order to study transverse and longitudinal BECs, we define variables in the LCMS ||19|| . 
Given a pair of particles, the LCMS is the frame of reference in which the sum of the two 
particle momenta, pi2 = (pi + P2), lies in the plane perpendicular to the thrust axis (see 
Fig. 0, where pu and p2i are the projections onto the plane). The momentum difference of 
the pair, Q = (p2 — Pi), computed in the LCMS, is resolved into the moduli of the transverse 
component, Qf, defined as shown in Fig. |^, and of the longitudinal component 

q^ = ipL-p;j^ (2) 

where the / direction coincides with the thrust axis. The momentum components are marked 
with a prime when they are measured in the LCMS. Qt may in turn be resolved into "out", 
Qtout^ and "side", Qt^e^ components 

Qt = Qwo + Qt,idc^ (3) 

where 6 and s are unit vectors in the plane perpendicular to the thrust direction, such that 
P12 = P120 defines the "out" direction and s = 1x6 defines the "side" direction. The symbol Q 
is used for the invariant modulus of the four-momentum difference Q = [IE2 — Ei|, (p2 — pi)]. 



It can be shown (see for instance ref. |2^) that, in the LCMS, the components Qi^.^^ and 
Q/ reflect only the difference in emission space of the two pions, while Qtout depends on the 
difference in emission time as well. Indeed, the scalar product between Q and the four-vector 
P = [(E2 + El), (p2 + pi)] is, in the LCMS 

Q ■ p = (e^-e;)(e'2 + e;) - q,„„,pi2. (4) 
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Since Q ■ P = 0, then 



Therefore, the BECs evaluated with respect to Qt^^^^ and in the LCMS yield information 
on the geometrical dimensions of the pion emitting source. In the string model LCMS 
represents the local rest frame of the string. In the following, we shall study the Bose-Einstein 
correlation function using two different definitions of the reference sample. 



4 The Bose-Einstein Correlation function 

The three-dimensional correlation function C is defined, in a small phase space volume around 
each triplet of Q,touti Qts»de Q' values, as the number of like-charge pairs in that volume 
divided by the number of unhke-charge pairs, used as a reference sample: 

^^TT+TT -'■^unlike 

Coulomb interactions between charged particles affect differently like- and unlike-charge 



pairs and thus modify the correlation function. A correction, based on the Gamow factors ||2T 
is applied: each pair of like-charge pions is weighted by a factor 

GKQ) = (e^^" - l)/27rr/ , (7) 

where rj = a^^m-^/Q,] each pair of unlike-charge pions is weighted by a factor 

G„(Q) = (1 - e-2-'')/2vrr/ . (8) 

The use of unlike-charge pairs as a reference sample gives large distortions of the correlation 
function in some regions of the domain (Qi„„t,Q4^.^^,Q;), caused by the presence of correlated 
n^'K~ pairs originating from decays of hadron resonances and of weakly decaying particles. 
In order to correct the correlation function, the Qt„„t, 0,1^^ and Q; distribution of the decay 
products of uj, rj, r]', K^, p", fo(980) and f2(1270) is determined for a sample of Jetset 7.4 [p^ 



multihadronic Monte Carlo events. Their contribution, due mainly to uj, rj, and p°, is 
subtracted from the unlike pion distribution NunUke- The OPAL version of Jetset used here 
reproduces the resonance structures reasonably well, although not perfectly. In particular the 
shape of the p° (0.64 < Q < 0.80 GeV) is not well modelled. For each resonance, differences 
between the simulated and measured |^ resonance rates are taken into account by scaling 
the distribution using the ratio of the measured production rate and the corresponding rate 
in Jetset. The resonance distributions obtained are then summed and scaled by the number 
of selected events. Figure |^ shows the distribution NunUke in the one-dimensional variable Q. 
The contribution of the decay products and the same distribution after the subtraction of this 
contribution are also shown. Two-dimensional projections of the function C(Qf„„j,Qi^.^^,Q;), 
after the correction for Coulomb and resonance decay effects, are shown in Fig. ^ for the sample 
of two-jet events selected, as an example, with jcut = 0.04 and with the third component, not 
plotted, limited to values up to 200 MeV. The histogram bin size, 40 MeV, is chosen to match 
the momentum resolution of the detector. The presence of Bose-Einstein correlations is seen 
as the sharp peak at small values of Qt„„„ Qt^ide Q'- 
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3500 




Q(GeV) 

Figure 2: The unlike-charge pairs distribution vs Q before (thick hne) and after (thin hne) the 
subtraction of the contribution from resonance decays (filled histograms; the higher peak is due 
to the the lower structure to the a;°, rj and K^j. 

5 Parametrisation of the Correlation function 



A minimum fit to the measured three-dimensional correlation function is performed using 



the following extended Goldhaber [24| parametrisation 



C(Qt„.„ Q*..,„ Qz) = N(l + Ae 



rf) 



(9) 



with 



The chaoticity parameter A measures the strength of the correlation, rt^„„ rf^.^^ and r; indicate 
the transverse and longitudinal extent of the two-pion source, N is a normalization factor 
necessary since the reference sample Nuniike is not normalized to the sample of like-charge 
pairs. The term F(Qi^^j, Q^^.^^, Q;) accounts for long-range two-particle correlations, due to 
energy and charge conservation and to phase space constraints. Alternative forms for this 
function have been considered, as in particular a fit with only the linear long-range terms 5jQj 
[i = iouti tside,^), and the results did not change significantly. The best results (lower values 
of x^/DoF and stability) are obtained with formula (10). The fits are performed over the 
range 0.04 < QtoutyQtsid^,Qi ^ 1-2 GeV. The region below 0.04 GeV is excluded because 
of the limited momentum resolution at low Qj values and of the presence of residual photon 
conversion pairs. Even after the subtraction procedure described in section 4, a few regions 
show significant distortions (see Fig. 3) assumed to be due to residual effects of pairs coming 



from resonance decays: these regions, corresponding to 0.28 < JQtout + Qts^de + Qz — 0-44 



GeV and 0.64 < JQt^J + Qt,,J + Qi < 0.80 GeV, are not used in the fits. 
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OPAL 




Figure 3: The projections of the three-dimensional correlation function C(Qt„^^,Qt^^^^,Qi), for 
two- jet events selected with jcut = 0.04, onto the (a) Qtout'Qi the (b) Qt^^^^-Qi planes when 
the third component has values up to 200 MeV. Notice the different behaviour of C as function 
of Qt^^a. a-nd Qi. 



Parameter 


Ycut = 0.01 


Ycut = 0.02 


Ycut = 0.04 


Ycut = 0.06 


inclusive sample 


N 


0.803 ±0.002 


0.815 ±0.001 


0.825 ±0.001 


0.831 ±0.001 


0.842 ±0.001 


A 


0.479 ±0.005 


0.464 ± 0.004 


0.454 ±0.004 


0.446 ± 0.004 


0.442 ±0.004 




0.522 ±0.007 


0.520 ±0.007 


0.525 ±0.006 


0.523 ±0.006 


0.536 ±0.006 


'^ts.d. [fm] 


0.750 ±0.007 


0.767 ±0.006 


0.783 ±0.006 


0.787 ±0.006 


0.809 ±0.006 


ri [fm] 


1.013 ±0.011 


1.014 ±0.010 


1.015 ±0.009 


1.011 ±0.009 


1.018 ±0.009 


K.. [GeV-i] 


-0.015 ±0.004 


-0.010 ±0.004 


-0.003 ±0.003 


-0.002 ±0.003 


0.004 ±0.003 


[GeV-i] 
.5/ [GeV-i] 


-0.149 ±0.004 


-0.148 ±0.004 


-0.146 ±0.003 


-0.145 ±0.003 


-0.142 ±0.003 


0.300 ±0.005 


0.262 ±0.004 


0.228 ±0.004 


0.212 ±0.004 


0.178 ±0.003 


[GeV-2] 


0.001 ±0.004 


0.010 ±0.003 


0.013 ±0.003 


0.015 ±0.003 


0.015 ±0.002 


[GeV-2] 
ei [GeV-2] 


-0.032 ± 0.004 


0.000 ±0.004 


0.022 ±0.003 


0.034 ±0.003 


0.053 ±0.003 


-0.072 ± 0.004 


-0.055 ±0.003 


-0.039 ±0.003 


-0.031 ±0.003 


-0.016 ±0.002 


xVDoF 


30089/24428 


31529/24428 


32758/24428 


33261/24428 


34632/24428 




1.351 ±0.027 


1.322 ±0.025 


1.296 ±0.021 


1.285 ±0.021 


1.258 ±0.020 



Table 1: Results of the fits of Eq. 9 to the measured three-dimensional correlation function 
C(QtoutMt,,de''^i) o^er the range 0.04 < Qt„„,,Qt,,^^,Q/ < 1.2 GeV, excluding the regions 
affected by residual resonance decay products described in the text. The quoted errors are 
statistical uncertainties obtained from the fits. The quality of the fits is indicated by the value 
ofx^DoF. 



The values of the parameters resulting from the fits to the correlation functions for two-jet 
events (with different values of jcut) and for the inclusive sample are given in Table |I|. For all 
the fits, the transverse and longitudinal radii are significantly different. In particular, the ratio 
"^ihtside between the longitudinal and transverse source radii is r^/r^^.^^ = 1.296 ± 0.021 [staf) 
for two-jet events selected with jcut = 0.04 and r^/r^^.^^ = 1.258 ± 0.020 {stat) for the inclusive 
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OPAL 




(a) (b) 

Figure 4: The projections of the three-dimensional correlation function C (Qt^^^,Qt^^j^,Qi), for 
two-jet events selected with jcut = 0.04, onto the (a) Qtout'Qi the (b) Qt^^^^-Qi planes when 
the third component has values up to 200 MeV. Notice the different behaviour ofC as function 
of Qu.,. and Qi. 



sample of events. 



6 The Correlation function C = C^^ta^qmc 

In order to improve the quality of the Goldhaber fits, reducing the effects from long-range 
correlations and resonance decay products, it is usual to define the following ratio of correlation 
functions, using data and Monte Carlo events: 



CDATA t\tDATA I t\tDATA 

nMC ~ AjMC lAjMC 
^ ^^like I unlike 



r^i ( r\ r\ \ Mfee / unlike !'\'\\ 

^ y'^touty^t.iae'^l) — r^Mr. ~ atMC /i\tMC • y^^) 



For this purpose, a sample of 7.2 million Jetset 7.4 multihadronic Monte Carlo events, which 
does not include BEC effects, is used. The Monte Carlo simulates many of the dynamical cor- 
relations present in the real data but not the Coulomb effect. Therefore, in Eq. 11, only N^^^-^"^ 
and ^ur^ike are corrected by the Gamow factors given in (7) and (8). The simulation includes 
the resonance decay products. The correction, due to the differences in the measured and sim- 
ulated resonance rates, is done by subtracting (adding) to the unlike-charge pairs sample in the 
Monte Carlo the fraction of pairs simulated in excess (deficit). Two- and one-dimensional pro- 
jections of the correlation function C, for the sample of two-jet events selected with jcut = 0.04, 
are shown in Fig. ^ and in Fig. ^, respectively. The Bose-Einstein correlations are clearly visible 
at small Qi„„t, Qt^idc Q'- The fits of Eq. 9, yield the parameters given in Table |[ As can be 
seen in the Table, the correlation function C is almost normalized; the slight difference from 1.0 
is due to the difference in the average multiplicity between data and MonteCarlo. The x^/DoF 
values for these fits are closer to unity than for those relative to the correlation function C. The 
dependences of rt^„„ rt^.^^, ri and of the ratios ri/rt^^^, rl/rt^^^^ on the jet resolution parameter 
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Figure 5: The projections of the three-dimensional correlation function C (Qt^^^,Qt^^^^,Qi), for 
two- jet events selected with jcut = 0.04, onto the (a) Qt^i^^ and (b) Qi axes when the two other 
components have values up to 200 MeV. 



Parameter 


Ycut = 0.01 


ycut = 0.02 


Ycut = 0.04 


ycut = 0.06 


inclusive sample 


N 


0.947 ±0.002 


0.950 ±0.002 


0.952 ±0.002 


0.954 ±0.001 


0.957 ±0.001 


A 


0.457 ±0.006 


0.446 ±0.005 


0.443 ±0.005 


0.441 ±0.005 


0.437 ±0.004 


rtaut M 


0.678 ±0.013 


0.654 ±0.012 


0.647 ±0.011 


0.637 ±0.010 


0.621 ±0.010 




0.781 ±0.010 


0.794 ±0.009 


0.809 ±0.009 


0.814 ±0.008 


0.831 ±0.008 


ri [fm] 


0.987 ±0.013 


0.989 ±0.012 


0.989 ±0.011 


0.989 ±0.010 


0.992 ±0.010 


^to„. [GeV-1] 


-0.004 ±0.005 


-0.002 ±0.004 


0.006 ±0.004 


0.005 ±0.004 


0.008 ± 0.003 


St..,. [GeV-i] 
Si [GeV-i] 


-0.062 ± 0.005 


-0.071 ±0.004 


-0.075 ±0.004 


-0.077 ±0.004 


-0.079 ±0.003 


0.114 ±0.005 


0.103 ±0.005 


0.092 ±0.004 


0.086 ±0.004 


0.073 ±0.003 


[GeV-2] 


0.012 ±0.004 


0.016 ±0.004 


0.012 ±0.003 


0.013 ±0.003 


0.011 ±0.003 


[GeV-2] 
ei [GeV-2] 


-0.020 ±0.005 


0.004 ±0.004 


0.019 ±0.004 


0.027 ±0.004 


0.037 ±0.003 


-0.045 ± 0.004 


-0.039 ±0.004 


-0.032 ±0.003 


-0.028 ±0.003 


-0.021 ±0.003 


xVDoF 


24654/24428 


25249/24428 


25398/24428 


25482/24428 


25836/24428 




1.264 ±0.033 


1.246 ±0.029 


1.222 ±0.027 


1.215 ±0.024 


1.194 ±0.024 



Table 2: Results of the fits of Eq. 9 to the measured three-dimensional correlation function 
C'(QtoutMt,,ae'Qi) over the range 0.04 < Qt<,„,,Qt,,^,,Qi < 1.2 GeV, excluding the regions 
affected by residual resonance decay products described in the text. The quoted errors are the 
statistical uncertainties obtained from the fits. The quality of the fits is indicated by the value 
ofx^/DoF. 



Ycut are shown in Fig. where the results for inclusive events are also presented. The main 
features of the results are a very slight decrease and increase, respectively, of r^^^^ and rt^.^^ 
as Ycut increases, while ri is independent of the jcut- The ratio ri/rt„„^ increases slightly, while 
ri/rt^^^^ decreases when jcut increases. The value of the chaoticity parameter A is between 0.457 
and 0.437 (decreasing slowly with jcut)- 
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Figure 6: Dependence on the resolution parameter jcut of (a) rt^^^, (h) rt^^^^ and (c) ri and of 
tie ratios (d) ri/rt^^^ and (e) n/rt^.^^ for the correlation function C (Qt„^^,Qt^.^^,Qi). 



The systematic uncertainties on the parameter values are estimated considering the de- 
viations with respect to a reference analysis, chosen to be the fit of Eq. 9 performed to the 
correlation function C {Qt^^^,Qt^.^^,Qi) for two-jet events selected with jcut — 0.04. In this case 
we have r^/r^^.^^ = 1.222 ± 0.027 (stat). The analysis is repeated changing some of the selection 
cuts: a maximum total momentum p < 30 GeV/c instead of p < 40 GeV/c and a charge 
unbalance smaller than 0.25 per event instead of 0.4. In order to check the stability of the 
results on the method used to correct the unlike-charge distribution for and resonance de- 
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cay products, the measured resonance rates are varied inside the experimental errors to obtain 
maximum and minimum values for the factors used to correct (i.e. either subtracted or added 
according to a defect or an excess with respect to the simulated resonance rates) N^^f-/^^. To 
estimate the influence of the long-range correlations, the fit is performed in a more restricted 
interval, 0.04 < Qf„„t, Qt^j^^,Qz < 1.0 GeV. Finally, the difference between the results of the 
fits to the correlation functions C and C is considered as an asymmetrical contribution to the 
systematic error. Table ^ shows the various contributions to the systematic error. The global 





n^ut [H 




n [fm] 




xVDoF 


a) Reference fit 


0.647 ±0.011 


0.809 ± 0.009 


0.989 ±0.011 


1.222 ±0.027 


25398/24428 


b) Modified track selection 

c) Max. resonance correction 

c) Min. resonance correction 

d) Restricted fit range 

e) Correlation function C 


0.656 ±0.012 
0.639 ±0.010 
0.657 ±0.011 
0.627 ±0.012 
0.525 ±0.006 


0.815 ±0.009 
0.812 ±0.008 
0.805 ± 0.009 
0.791 ± 0.009 
0.783 ± 0.006 


0.995 ±0.012 
0.988 ±0.010 
0.989 ±0.011 
0.975 ±0.011 
1.015 ±0.009 


1.221 ±0.028 
1.217 ±0.024 
1.228 ±0.027 
1.233 ±0.028 
1.296 ±0.021 


25209/24428 
25566/24428 
25319/24428 
13751/13053 
32758/24428 



Table 3: Results of the fit of Eq. 9 to several variations of C (Qtout^Qts^de^ Q-O' ^ listed in the 
text and in the first column. The quoted errors are only statistical. 



systematic uncertainties are computed by adding in quadrature the differences between the 
reference fit a) and the variations b) - e). 

The conclusion from this analysis is that, as observed in the LCMS, the pion emitting region 
is elongated, with r^ greater than r^^^^^. From Fig. 6 it is evident that the ratio ii/it^^^^ has a 
(small) dependence on jcuu the largest value is obtained for smaller jcut- One also observes 
that r^ is independent on jcut, while r^^.^^ increases with increasing jcut- As an example, we 
quote the following parameter values obtained for two-jet events, selected using jcut = 0.04: 

ri„„, = (0.647 ± 0.011 (stat) t^^H (syst)) fm 
Tt^^,^ = (0.809 ± 0.009 (stat) 1°;°^^ (syst)) fm 
r, = (0.989 ± 0.011 (stat) (syst)) fm 

(12) 

n/it^^^^ = 1.222 ± 0.027 (stat) t";"!^ (syst). 



The results of this analysis are in qualitative agreement with recent results from the L3 
collaboration |jl5[. In the L3 analysis, which uses an event-mixing technique to compute the 
reference sample, the ratio of transverse to longitudinal radius is 0.81 ± 0.02 toil , correspond- 
ing to r,/r,^^,^ = 1.23 ± 0.03 t'^f, . 

The results can also be compared with the predictions of a recent model of BECs based 
on string fragmentation |T^. In this model, the different mechanisms that generate the longi- 



tudinal (i.e. along the string) and transverse momenta of the particle, lead to a longitudinal 
correlation length, representing the space-time difference along the string between the produc- 
tion points, which is larger than the transverse correlation length. 
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Figure 7: The one- dimensional correlation function C'(Q) with the result of the one- dimensional 
Goldhaber fit, Eq. 13, superimposed. 



As a check, a two-dimensional analysis where the transverse component Qt of the pair mo- 
mentum difference is not split into "out" and "side" components, is done. The two-dimensional 
Goldhaber function fitted to the correlation function C'(Qt,Qi) gives, in the case of two-jet 
events selected with jcut = 0.04, r^/r^ = 1.360 ± 0.026 (stat), in agreement with a longitudinal 
source size larger than the transverse size. While r^^^^ and Tt^.^^ show a slight dependence on 
the jet resolution parameter jcut (in opposite directions, see Fig. 6), the parameter it obtained 
from the two-dimensional analysis is independent on jcut- 

We conclude that we always have ii/it^e greater than one. 

7 One-dimensional analysis of the inclusive sample 

The one-dimensional correlation functions C(Q) and C'(Q), where Q is the modulus of the 
four-momentum difference of the pion pair, are studied for the inclusive sample. The results 
can be compared with those published in p|; the present analysis uses more data (4.3 million 
TP hadronic decays instead of 3.6 millions) and a different version of the Jetset Monte Carlo 
(7.4 instead of 7.3). The correlation function C'(Q), after the corrections for Coulomb and 
resonance decay products, is shown in Fig. |^. The one-dimensional Goldhaber function 

C'(Q) = N{1 + Ae-Q''"') (1 + 5Q + eQ^) (13) 

is fitted to the measured correlation function over the range 0.04 < Q < 1.2 GeV. There 
are apparent distortions in the unlike-charge pairs distribution even after the correction for 
the contribution of resonance decays products, as discussed in the three-dimensional analysis. 
Therefore the Q regions affected most by tu, ?7, K^., p° and /o(980) decay products, corresponding 
respectively to 0.28 < Q < 0.48 GeV, 0.68 < Q < 0.84 GeV and 0.96 < Q < 1.0 GeV, are 
excluded from the fit range. 
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r [fm] 


A 


xVdoF 


a) Reference fit 


i.002 ± 0.0i6 


0.574 ± 0.0i9 


69/15=4.6 


b) Modified track selection 

c) Max. resonance correction 

c) Min. resonance correction 

d) Restricted fit range 

e) Correlation function C 


0.984 ± 0.0i7 
0.99i ± 0.0i6 
i.0i4 ± 0.0i7 
l.Olf ± 0.023 
0.909 ± O.Oi? 


0.590 ± 0.02i 
0.605 ± 0.019 
0.546 ± 0.019 
0.557 ± 0.024 
0.580 ± 0.022 


55/15=3.7 
66/15=4.4 
75/15=5.0 
67/10=6.7 
150/15=10.0 



Table 4: Results of the fits of Eq. 13 to several variations of C'(Q). The quoted errors are only 
statistical. 



The systematic error on the measured values of the parameters r and A is evaluated in 
the same way as done for the three-dimensional analysis. The results of the various fits are 
summarized in Table |[ 

The fit gives the following values for the parameters: 

r = (1.002 ± 0.016 (stat) toml (syst)) fm , A = 0.574 ± 0.019 (stat) 1°;°^^ (syst). (14) 

Since the percentage number of charged tracks which are pions is about 87%, one can esti- 
mate that the value of the A parameter would be a factor of 1.32 larger in the case of a 100% 
pure pion sample. 

The values of Eq. 14 are in agreement and replace the values previously published by the 
OPAL Collaboration, see ref. [@]. 



8 Conclusions 

Using 4.3 million hadronic events from Z° decays, the Bose- Einstein correlation function for 
two identical charged bosons, mainly vr^vr^, is studied in the three components of the mo- 
mentum difference, longitudinal and transverse ("out" and "side" components) with respect 
to the thrust direction, in the Longitudinally CoMoving System. The geometrical structure 
of the source is obtained from an extended Goldhaber fit of Eq. 9 to the Coulomb corrected 
BEC functions C{Qtout:Qt,,aeyQi) and C'(Qt,„t,Qt,,d,,Qz). In all cases the longitudinal radius is 
significantly larger than the transverse radius. 



The longitudinal and transverse radii of the emitting region are studied as a function of the 
two-jet resolution parameter jcut- The analyses indicate that, as ycut increases, r^^.^^ increases 
slowly, r; remains constant and that the ratio ii/itMe decreases. In the framework of the string 
model of ref. |jl3|, the observed different values of the transverse and longitudinal correlation 
lengths are explained in terms of two different generation mechanisms of the longitudinal and 
transverse momentum components with respect to the string direction. 



The fit of Eq. 9 to C'(Qt^^^,Q(^.^^,Q;), in the case of two-jet events selected with jcut = 0.04, 
yields the parameter values given in (12), in particular 11/11^^ ~ 1-222 ± 0.027 [stat) l^oii 
{syst). The corresponding value for the inclusive sample is ii/itMe. ~ 1-194, with similar uncer- 
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tainties. 



The inclusive sample of events is also analysed in terms of the one-dimensional correlation 
function C'(Q). The fit gives the parameter values r = (1.002 ± 0.016 (stat) to^m i^y^i)) 
and A = 0.574 ± 0.019 (stat) {syst). 

In conclusion, the present analysis shows that the emitting source of two identical pions, 
measured in e^e~ interactions at energies close to the Z° peak, is elongated. In particular, as 
computed in the LCMS, the emitting source for two identical pions has global dimensions of 
about 1 fm, but with the longitudinal dimension about 20% larger than the transverse dimen- 
sion. 
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